-----------------------------------------------------------------------------------------------------------------------------------
Get this COARE soundings .nc data file It is 1726 soundings on a 5mb grid 1000-100 mb. (181 values)
Get the IDL or Matlab codes: main ones are esat, mixrat (or saturation mixing ratio), entropy, invert_entropy
Read the COARE soundings into memory. Notice array sizes.
%------------------------------------------------------------ % get data : array [181,1726] %------------------------------------------------------------ p = nc_varget (fname, 'p'); % P is constant 5mb interval, 1000-100.[1000,995,990,...110,105,100] T = nc_varget (fname, 'T'); RH = nc_varget (fname, 'RH'); % [181,1726] irtemp = nc_varget (fname, 'IRTEMP'); % [1726] %------------------------------------------------------------------
Compute q (mixing ratio) = {[0.622 *esat(T) / (p-esat(T)) ] *RH} /100 = mixingratio(esat(T),p) *RH / 100
% #4. compute q : traditionally w and ws are used for q (q is the same as qv) and qsat % ws = 0.622 .* es ./ (press-es); % w=(ws.*rh)/100.; % what this value means? kg/kg
*Figures : mean profile of T, q(and q_sat) and RH
Compute the entropy, total water, and pressure of a near-surface parcel. Call these s0, qt0, p0. %=====================================================================
% #6. Compute s0, qt0, and p0
%---------------------------------------------------------------------
T0=T_mean(1);
p0=p_mean(1);
qt0=q_mean(1);
s0=entropy(T0, p0, qt0);
% entropy of the actual sounding (qtot=q=qv, b/c no condensate)
ent = entropy(T,p,q);
% SATURATION ent of the sounding (qtot=qsat)
entsat = entropy(T,p,mixingratio(esat(T),p));
% dry version (qtot=> q=0)
entdry = entropy(T,p,1e-7*q);
Figures: mean profiles of Entropy
Same as previous figure, but y axis is height in this figure.
matlab code
-----------------------------------------------------------------------------------------------------------------------------------
Get this COARE soundings .nc data file It is 1726 soundings on a 5mb grid 1000-100 mb. (181 values)
Get the IDL or Matlab codes: main ones are esat, mixrat (or saturation mixing ratio), entropy, invert_entropy
Read the COARE soundings into memory. Notice array sizes.
%------------------------------------------------------------
% get data : array [181,1726]
%------------------------------------------------------------
p = nc_varget (fname, 'p'); % P is constant 5mb interval, 1000-100.[1000,995,990,...110,105,100]
T = nc_varget (fname, 'T');
RH = nc_varget (fname, 'RH'); % [181,1726]
irtemp = nc_varget (fname, 'IRTEMP'); % [1726]
%------------------------------------------------------------------
Compute q (mixing ratio) = {[0.622 *esat(T) / (p-esat(T)) ] *RH} /100 = mixingratio(esat(T),p) *RH / 100
% #4. compute q : traditionally w and ws are used for q (q is the same as qv) and qsat
% ws = 0.622 .* es ./ (press-es);
% w=(ws.*rh)/100.; % what this value means? kg/kg
*Figures : mean profile of T, q(and q_sat) and RH
Compute the entropy, total water, and pressure of a near-surface parcel. Call these s0, qt0, p0.
%=====================================================================
% #6. Compute s0, qt0, and p0
%---------------------------------------------------------------------
T0=T_mean(1);
p0=p_mean(1);
qt0=q_mean(1);
s0=entropy(T0, p0, qt0);
% entropy of the actual sounding (qtot=q=qv, b/c no condensate)
ent = entropy(T,p,q);
% SATURATION ent of the sounding (qtot=qsat)
entsat = entropy(T,p,mixingratio(esat(T),p));
% dry version (qtot=> q=0)
entdry = entropy(T,p,1e-7*q);
Figures: mean profiles of Entropy
Same as previous figure, but y axis is height in this figure.
working on #12.